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Abstract. We show that the sharp Young's inequality for convolutions first 
obtained by Bechner [2] and Brascamp-Lieb [7] can be derived from the mono- 
tone in time evolution of a Lyapunov functional of the convolution of two 
solutions to the heat equation, with different diffusion coefficients, first intro- 
duced by Bennett and Bez in [4]. Our proof is based on a suitable adaptation 
of an old idea of Stam [19] and Blachman [6], used to obtain Shannon's 
entropy power inequality. 
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1 Introduction 

The goal of this note is to present a new proof of the Young inequality in the sharp 
form obtained by Bechner [2], 



\\f*g\\ r <{A p A q A r ,) n \\f\\ p \\g\ 



q- 



In (1) / E L p (R n ), g E L q (R n ), 1 < p, q, r < oo and 1/p + 1/q = 1 + 1/r. Moreover, 
the constant A m which defines the sharp constant is given by 

/ m l/m \ V2 

A ™ = ( _,i/ m / ) ( 2 ) 



where primes always denote dual exponents, l/m + l/m'=l. 

The best constants in Young's inequality were found by Beckner [2], using ten- 
sorisation arguments and rearrangements of functions. In [7], Brascamp and Lieb 
derived them from a more general inequality, which is nowadays known as the 
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Brascamp-Lieb inequality. The expression of the best constant, in the case in which 
both / and g are probability density functions, is obtained by noticing that inequal- 
ity (1) is saturated by Gaussian densities. This principle has been largely utilized 
by Lieb in a more recent paper [15]. Among many other results, this paper contains 
a new proof of the Brascamp-Lieb inequality. In [7], Brascamp and Lieb noticed 
that the sharp form of Young inequality also holds in the so-called reverse case 

\\f * g\\ r > (A p A q A rl ) n \\f\\ p \\g\U, (3) 

where now < p, q, r < 1 while, as in Young inequality (1), 1/p+l/q — 1 + 1/r. In 
this case, however, the dual exponents p', q', r' are negative, and 

/ »"•' \ 1/2 

The proof of this sharp reverse Young inequality was subsequently simplified by 
Barthe [1]. While the original proof in [7] was rather complicated, and used ten- 
sorisation, Schwarz symmetrization, Brunn-Minkowski and some not so intuitive 
phenomenon for the measure in high dimension, the new proof in [1] was based on 
relatively more elementary arguments and gave a unified treatment of both cases, 
the Young inequality (1) and its reverse form (3). As a matter of fact, the proof of 
the main result in [1] relies on a parametrization of functions which was used in [13] 
and was suggested by Brunn's proof of the Brunn-Minkowski inequality. 

In a recent paper, Young's inequality has been seen in a different light by Ben- 
nett and Bez [4]. In their paper, Young's inequality is derived by looking at the 
monotonicity of a suitable functional of the solution to the heat equation. Even if 
not explicitly mentioned in the paper, this idea connects Young's inequality in sharp 
form with other inequalities, for which the proof exactly moved along the same idea. 

The connections of the sharp form of Young inequality with other inequalities 
has been enlightened by Lieb in [14]. He proved in fact that, by letting p, q, r — > 1 in 
(1), the sharp form of Young's inequality reduces to another well-known inequality 
in information theory, known as Shannon's entropy power inequality [18]. 

In its original version, Shannon's entropy power inequality gives a lower bound 
on Shannon's entropy functional of the sum of independent random variables X, Y 
with densities 

exp (2H(X + Y)) > exp (2H(X)) + exp (2H(Y)) , (5) 

with equality if X and Y are Gaussian random variables. Shannon's entropy func- 
tional of the probability density function f(x) of X is defined by 

H{X) = H(f) = - f f(x) log f(x) dx. (6) 
Jr 
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Note that Shannon's entropy functional coincides to Boltzmann's if-functional [10] 
up to a change of sign. The entropy-power 

N(X)=N(f)=exp(2H(X)) 

(variance of a Gaussian random variable with the same Shannon's entropy func- 
tional) is maximum and equal to the variance when the random variable is Gaussian, 
and thus, the essence of (5) is that the sum of independent random variables tends 
to be more Gaussian than one or both of the individual components. 

The first rigorous proof of inequality (5) was given by Stam [19] (see also Blach- 
man [6] for the generalization to n-dimensional random vectors), and was based on 
an identity which couples Fisher's information with Shannon's entropy functional 
[12]. 

The original proofs of Blachman and Stam make a substantial use of the solution 
to the heat equation 

^ = A/(:M), (7) 

that is, for t > 0, of the function f(x,t) = f * M 2 t(x), where M t (x) denotes the 
Gaussian density in W 1 of variance t 

M ^) = 7tA^^p( 1 4-)- (8) 



(27rt)»/ 2 ' V 2* 

Other variations of the entropy-power inequality are present in the literature. Costa's 
strengthened entropy-power inequality [11], in which one of the variables is Gaus- 
sian, and a generalized inequality for linear transforms of a random vector due to 
Zamir and Feder [22]. 

Also, other properties of Shannon's entropy-power N(f) have been investigated 
so far. In particular, the concavity of entropy power theorem, which asserts that 

^ 2 (N(f*M 2t ))<0 (9) 

Inequality (9) is due to Costa [11]. More recently, a short and simple proof of (9) 
has been obtained by Villani [21], using an old idea by McKean [17]. 

Summarizing, the proof of Stam is based on the following argument. Let f(x,t) 
and g(x, t) be two solutions of the heat equation (7) corresponding to the initial data 
f(x) (respectively g(x)). If the entropies of the initial data are finite, one considers 
the evolution in time of the functional Qf, g (t) defined by 

Evaluating the time derivative of Qf t9 (t), and using a key inequality for Fisher 
information on convolutions, shows that 0/, s (t) is increasing in time, and converges 
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towards the constant value 0/ )9 (+oo) = 1, thus proving inequality (5). Note that 
this method of proof also determines the cases of equality in (5). It is interesting to 
remark that the evaluation of the limit of 0/ iS (t), as t — > oo, is made easy in reason 
of a scaling property. Indeed, the (Lyapunov) functional Q(f,g) is invariant with 
respect to the scaling 

= (Jr) ' a>0 ' ( n ) 

which preserves the total mass of the function /. The importance of this property 
will be clarified later on. 

The proof by Stam is a physical proof, in the spirit of Boltzmann iJ-theorem 
[10] in kinetic theory of rarefied gases, where convergence towards the Maxwellian 
equilibrium is shown in consequence of the monotonicity in time of the logarithmic 
entropy (6). 

In the rest of the paper, inspired by the Stam's approach to the proof of Shan- 
non's entropy power inequality, we will present a physical proof of both direct and 
reverse Young's inequalities, which is based on the two ingredients specified above: 
a suitable use of solutions to the heat equations, coupled with the scaling invariance 
property (11). Our proof is alternative to the proof of [4], and relies on a result 
which generalizes Stam proof of subadditivity of Fisher information. 

In the following Section 2 we will describe how the method works by proving 
Holder inequality. Even if this is a well-known result [8, 4, 3], it will give indication 
on the underlying methodology. Next, Section 3 will contain our main result. 

To make computations as simple as possible, we will present all proofs in one 
dimension. Without loss of generality, in fact, one can easily argue that identical 
proofs hold in dimension n, with n > 1. Also, if not strictly necessary, we will 
restrict ourselves to consider as functions probability density functions. Since our 
computations involve solutions to the heat equation, all details involving regularity 
and the various integrations by parts can be assumed to hold true. 

The basic idea used here is that many inequalities can be viewed as the tendency 
of various Lyapunov functionals of the solution to the heat equation to reach their 
maximum value as time tends to infinity. The discovery of a Lyapunov functional 
which allows to prove Young inequality is only one of the possible application of this 
idea [3, 4, 5]. In a forthcoming paper [20], we are going to develop this strategy by 
revisiting various well-known inequalities in terms of the monotonicity of suitable 
Lyapunov functionals. 

2 Heat equation and Holder inequality 

We begin by showing that Holder's inequality can be viewed as a consequence of 
the time monotonicity of a suitable Lyapunov functional of the solution to the heat 
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equation [8, 4]. 

Holder's inequality for integrals states that, if p, q > 1 are such that 1/p+l/q = 1 

/ \f(x)g(x)\dx<( [ \f(x)\*dx) ^ ([ \g(x)\*dx) '* . (12) 

Moreover, there is equality in (12) if and only / and g are such that there exist 
positive real numbers a and b such that af p (x) = bg q (x) almost everywhere. Holder's 
inequality can be proven in many ways, for example resorting to Young's inequality 
for constants, which states that, if 1/p +l/q — l 

C d < - + _ (13) 
V Q 

for all nonnegative c and d, where equality is achieved if and only if c p = d q . 

Without loss of generality, one can assume that the functions /, g in (12) are 
nonnegative. A different way to achieve inequality (12) is contained into the follow- 
ing 

Theorem 1. Let $ u>v (t) be the functional 

*„,„(*)= / u(x,ty'*v(x,t) 1/q dx, (14) 

where 1/p+l/q = 1, andu(x,t) and v(x,t) , t > 0, are solutions to the heat equation 
corresponding to the initial values u(x) G L : (R) (respectively v(x) G L 1 (R)j. Then 
<f>u,v(t) is increasing in time from 



= 0)= / uixf^v^x) 1 ^ dx, 



to 

lim & u ,v(t) = I u(x) dx I v(x)dx. 



t— >OD 



Proof. We outline that the functional $ u ,v(t) is invariant with respect to the scaling 
(11). Moreover, the condition u(x),v(x) G L 1 (IR) is enough to ensure that & UiV (t) G 
L 1 (M) at any time t > 0. Indeed, inequality (13) implies 



u(x,t) l/p v(x,t) llq < -u(x,t) + -v(x,t), 
p q 



where, since u(x, t) and v(x,t) are solution to the heat equation, 

/ u{x ) t)dx= I u(x)dx, / v(x,t)dx= / v(x) dx. 
Jm. Jr. Jr Jr 



Let us first proceed in a formal way. However, by resorting to the smoothness 
properties of the solution to the heat equation, all mathematical details can be 
rigorously justified. 

Let us evaluate the time derivative of $(£). It holds 

*'«,«(*)= / [WM) 1/p )MM) 1/? + u(M) 1/p WM) 1/? )J dx = 
Jr 

[ -u l/p - 1 v 1/q u xx + -u 1/p v l/q - 1 v xx dx= I -u~ 1 / q v l/q u xx + -u 1/p v- 1/p v 
JrLP q J Jr LP q 

Integrating by parts we end up with 



dx. 



pq Jr 



^l) 2 - 2—— + (—\ 
u J U V V V ) 



dx 



1 

Pq Jr 



u 



i/p, 



x,t)v 1/q (x,t) 



u x (x,t) v x (x,t) 



dx > 0. 



(15) 



u(x,t) v(x,t) 

Hence the functional ® UiV (t) is increasing in time. Note that the time derivative of 
the functional is equal to zero if and only if, for every t > 

u x (x,t) _ v x (x,t) = 
u(x,t) v(x,t) 

for all points iGl This condition can be rewritten as 



d j u(x,t) 
dx v(v,t) 



0. 



Consequently $'(t) = if and only if 

u(x, t) = cv(x, t) 



(16) 



for some positive constant c. Thus, unless condition (16) is verified almost every- 
where at time t — 0, the functional $(t) is monotone increasing, and it will reach 
its eventual maximum value as time t — > oo. The computation of the limit value 
uses in a substantial way the scaling invariance of $. In fact, at each time t > 0, 
the value of & u>v (t) does not change if we scale u(x,t) and v(x,t) according to 



u(x, t) ->■ U(x,t) = Vl + 2tu(xVl + 2t,t) 
v(x,t) ->■ V(x,t) = Vl + 2tv(xVl + 2t,t). 

On the other hand, it is well-known that [9] 



(17) 



lim U(x,t) = Mi (a;) u(x) dx ]imV(x,t) = M 1 (x) v(x)dx, (18) 



where, according to (8) Mi(x) is the Gaussian density in R of variance equal to 1. 
Therefore, passing to the limit one obtains 

Urn = (^J u{x)dx^j (^J v(x)dx^j J M 1 (x) 1/p M 1 (x) 1/q dx = 

u(x)dx) I / v(x)dx) J Mi(x) dx — ( / u(x)dx) ( / dx ) 

Since 



lim $ uv (t) = [ u{x) 1/p v(x) 1/q dx, 
the monotonicity of the functional $(t) implies the inequality 

/ u(x) 1/p v(x) 1/q dx < ( [ u(x)dx) 1 ( [ v(x)dx) '\ (19) 

with equality if and only if (16) is verified at time t — 0, that is 

u(x) = cv(x), (20) 

for some positive constant c. Setting / = u 1 ^ and g = v x l q proves both Holder 
inequality (12) and the equality cases. □ 

Despite its apparent complexity, this way of proof is based on a solid physical 
argument, namely the monotonicity in time of a Lyapunov functional of the solution 
to the heat equation. This gives a clear indication that many inequalities reflect 
the physical principle of the tendency of a system to move towards the state of 
maximum entropy. In the next Section we will see how this idea applies to prove 
Young's inequality. 



3 Young's inequality and Lyapunov functionals 

The proof of the sharp Young's inequality follows along the same lines of the proof 
of Holder's inequality we presented in Section 2. In this case the key functional to 
study is the one considered by Bennett and Bez [4] 

*«,«(*) = (j[ (MM) 1/p *^(M) lA? ) r dx^J ' , (21) 

where, as in Young's inequality, l/p+ 1/q — 1 + 1/r. With respect to the notations of 
the previous Section, there is a substantial difference in the meaning of the functions 
u(x,t) and v(x,t). Here u(x, t) and v(x,t) are still solutions of the heat equation 
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corresponding to the initial data u(x) (respectively v(x)). However, these solutions 
correspond to two different heat equations, with different coefficients of diffusions, 
say a and j3. In other words, u(x,t) solves the diffusion equation 

ut = ®Uxx, (22) 

while v(x, t) solves 

v t = f3v xx . (23) 

Hence u(x,t) and v(x,t) diffuse at different velocities. It is a simple exercise to 
verify that, in view of the relationship between p,q and r, the functional *& UtV (t) is 
invariant with respect to the mass preserving scaling (If). 

Theorem 2. Let ^ UjV (t) be the functional (21), where l/p+ l/q = 1 + 1/r, and 
u(x,t) and v(x,t), t > 0, are solutions to the heat equation corresponding to the 
initial values u(x) G L 1 (M) (respectively v(x) G L 1 (M,)). Then, ifp,q,r > 1, and the 
diffusion coefficients in (22) and (23) are given by a = q'/p (respectively (3 = p'/q), 
or by a multiple of them, ^ UjV (t) is increasing in time from 



1/r 



to the limit value 

Urn V v , v (t) = (A p A q A r ,) 1/2 (^Ju(x)dx S j (^Jv{x)dx S j . (24) 



// on the contrary < p, q, r < 1, and the diffusion coefficients in (22) and (23) 
are given by a = \q'\/p (respectively (3 = \p'\/q), or by a multiple of them, ^ u , v (t) is 
decreasing in time from 



= 0) = ^ (u(x) 1 / p *v(x) 1 ' q Y dx 



1/r 



to the limit value (24), where now A m is given by (4)- In both cases W uv (t) = if 
and only ifu(x,t) andv(x,t) are Gaussian. 

Proof. Let us consider first the case in which p,q,r > 1. Without loss of generality, 
we will assume that both the initial data u(x) and v(x) are probability density 
functions. This is sufficient to show that, for any time t > 0, the functional ^ u ,v{t) 
is bounded. Indeed we can write 



/ 



u{x-y) llp v{y) l/q dy 
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/ u(x - y) 1/p v{y) 1/q dy + I u(y) 1/p v(x - y) 1/q dy. 

J {u(x-y)<v(y)} J {u(y)>v(x~y)} 

Now, since r > 1, and v(x, t) has mass equal to 1, Jensen's inequality implies 

(I u(x - y) 1/p v{y)^ q dy) = 

\J \u(x-y)<v(y)} J 



J 

J {u(x— 



{u(x-y)<v(y)} 

uix-yf/vviyfl^v^dy) < 



y)<u(y)} 



Note that 



/ u(x - y ) r/p v(y) r/q ~ r v(y) dy 

J {u(x-y)<v(y)} 

I u(x-y) r/p v(y) r/q - r+1 dy. 

J {u(x-y)<v(y)} 



rj-t rj-t rjri 

- + -- r + l = 2, - > 1. 
p q p 



Therefore, on the set {u(x — y) < v(y)}, since the exponent of v (y) is smaller than 
1, 

u(x - y) r/p v(y) r/q - r+1 < u(x - y)v{y). (25) 
Inequality (25) follows simply dividing by u 2 . Therefore 

I u ( x _ yy/Pviy) 1 ^ dy) < j u{x-y)v{y)dy<\. 

' {u(x-y)<v(y)} J J {u(x~y)<v(y)} 

Identical computations show that 

(I u{y) l/p v{x - yf/ q dy) < [ u(y)v(x - y) dy < 1. 

\J {u(y)>v(x-y)} J J {u(x-y)<v(y)} 

Therefore 

{u(x) l ' p *v(x) l ' q ) r < 2c r , 



where c r is the positive constant in the inequality 

(a + b) r < c r (a r + b r ). 

We proceed now to compute the time derivative of the functional ty UtV (t). To shorten, 
let us denote 

h(x, t) = u(x, t) 1/p * v(x, t) l/q . (26) 
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Since u(x,t) and v(x,t) are solutions to the heat equation 

^h(x, t) = j t j^u(x- y, t)^v(y, t) 1 '* dy = 

- [ u (x-y,t) 1 ^ 1 u t (x-y,t)v(y,t) 1 ^dy + - [ u{x-y,t)^v{y,tf^ l v t {y,t) dy = 

pJr q Jr 

[ u(x-y,t) 1 / p - 1 u xx (x-y,t)v(y,t) 1 ^dy+^- [ uix-y^fl^tf^v^t) dy. 
Jr Q Jr 



V 
a 

V 

On the other hand, we have 
d 2 



dx 2 



Hence 



h(x, t) = J^ (^u(x - y, tf/^u^x - y, t)j v(y, t) 1 '* dy = 

- [ u(x-y,t) 1 / p - 1 u xx (x-y,t)v(y,t) 1 /< 1 dy+ 

V Jr 

-f--l) [ u(x- y ,t)^- 2 u 2 x (x- y ,t)v( y ,t)^d y . (27) 
V \V J Jr 

- / u(x-y,t) 1/p - 1 u xx (x-y,t)v(y,t) 1/q dy = 

V Jr 

h(x,t) + ± [ u(x - y ,t)^v(y,t) 1/q (-Y (x -y,t)dy. (28) 
VP Jr v u J 



dx 2 

Analogous formula for the last integral in (27). Therefore we have 



3 d 2 
Qjh{x,t) = (a + p)—h(x,t)+ 



a 



. u{x - y, t)^v(y, ty/o (^) 2 (x - y, t) dy+ 
VV Jr v u J 

j u(x - y, ty/ p v(y, t) 1/q (^) 2 (v, t) dy. (29) 



JL 

qq' 

Making use of formula (29), we obtain 



d 
di 



h r (x,t)dx = r / h(x,t) r 1 h t (x : t) dx = r(a + (5) / h r 1 (x, t)h xx (x, t) dx+ 
Jr Jr 

[ h r ~\x,t) \— t f u(x-y,ty/vv(y,t) 1/q (-y (x-y,t)dy+ 
Jr IVV Jr \ u / 

A / u(x-y,t)^v(y,t)^ ( V ^)\y,t)dy 
QQ Jr \ v / 
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Since it holds 

y}/P p u ' v Vp q V 

we obtain 

2 



and 

ft f /?) \2 « /• (^^(y^)) 

7 / "(s-v^'Mv,*) 1 ' 9 ^) G/,*) = / n^-y) K vVq J d y . 

Jr yv/ 9 Ji v 1 yy) 



qq 
Finally 



/ h r (x,t)dx = -(a + p)(r- 1) / /i r - 2 (x,t) (h x (x,t)) 2 dx+ 

a L [ h r -\x,t)A(u 1/p ,v 1/q )(x,t)dx + (3^- [ h r -\x,t)B(u 1 / p ,v 1/q )(x,t)dx. (31) 
P Jr q Jr 

In (31) we denned 

A(f,g)(x,t) = [ ifx( f- y ^ )2 g(y,t)dy, (32) 
Jr j\ x ~y) 

and 

W, </)(*, f) = jf /(* - dy. (33) 



Since 



dt 

the sign of the time derivative of the functional ^ u ,v(t) depends of the sign of the 
expression on the right-hand side of (31). In order to determine this sign, the 
following Lemma will be of paramount importance. 

Lemma 3. Let f(x) and g(x) be probability density functions such that both A(f,g) 
and B(f, g), given by (32) and (33), are well defined. Then, for all positive constants 
a, b and r > 

(a 2 + b 2 + 2abr) [ (f * gf' 2 ((/ * g) x ) 2 dx < 
Jr 

a 2 f(f* gy-'Aif, g) dx + b 2 I \f * gf^Bif, g) dx. (34) 
Jr Jr 

Moreover, there is equality in (34) if and only if, for any positive constant c and 

constants m\,m<i, f and g are Gaussian densities, f(x) = M ca [x — mi) and g(x) = 

M cb (x - ra 2 ). 
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Proof. The proof will follow along the same lines of the analogous one for Fisher 
information, given by Blachman [6]. First of all, to easily justify computations, let 
us prove the lemma by considering smooth functions / and g. Then the proof for 
general / and g will follow owing to the convexity properties of A and B [16]. This 
can be easily done by considering f * M t and g * M t solutions to the heat equation 
for some t > 0. Let 

k(x) = f*g(x). 
Then, for any pair of positive constants a, b 

(a + b)k'(x) = a / f'(x-y)g(y)dy + b / f(x-y)g'(y)dy. 
Jr Jr 



Therefore 



(a + b) 



k'(x) 
k(x) 



f'{x-y) f(x-y)g(y) 
f(x-y) k(x) 



dy + b 



g'iv) f{x-y)g{y) 

g{y) k(x) 



dy = 



where we denoted 



[ ( f'(x- y) , J(v)\ . , , 



Note that, for every x G M, d\i x is a unit measure on R. Consequently, by Jensen's 
inequality 



(a + bf 



k\x) 
k(x) 



f'(x-y) g'(y) 
a— r + b 



f(x-y) g(y) 



d^xiy) 



< 



(35) 



Hence, for every constant r > 

(a + bf I 

Jm. 

k r (x 
k r ~\x) 



k r (x) 



k'(x) 
k(x) 



dx < 



f'(x-y) ,g'(y)\ 2 f(x-y)g(y) 



[ k r ~\x) a 2 I 



f(x-y) +b g(y) 

U\x-y)f 



f(x-y) 



g(y) dy + b' 



k(x) 

(g'(y)) 2 
g(y) 



dy dx = 



f{x-y) dy 



dx+ 



2ab [ k r ~\x) [ f'(x-y)g'(y)dydx. 
Jr Jr 
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On the other hand, 

/ f'(x-y)g'(y)dy = k"(x), 
Jr 

so that 

/ k r ~\x) [ f'(x-y)g'(y)dydx = 
Jr Jr 



k r (x)r"(x) dx = -(r - 1) / k r -\x)(k\x)) 2 dx. 

Jr 

This concludes the proof of the lemma. The cases of equality are easily found from 
the following argument. Equality follows if, after application of Jensen's inequality, 
there is equality in (35). On the other hand, for any convex function tp and unit 
measure djj, on the set Q, equality in Jensen's inequality 

<p(f fdfM)< f <p(f)dpi 
Jn Jn 

holds true if and only if / is constant, so that 

/= / fd/M. 
Jn 

In our case, this means that there is equality if and only if the function 

a f'(*-v) | h 9\y) 

f(x-y) g(y) 

does not depend on y. If this is the case, taking the derivative with respect to y, 
and using the identity 

d ff{x-y)\ d ff'(x-y) 



dy V f( x ~y) ) dx \f(x -y) J ' 
we conclude that / and g have to satisfy 

d 2 d 2 
a dx^ l ° gf ( X ~^ = b df l0g9 ^' (36) 

Note that (36) can be verified if and only if the functions on both sides are constant. 
Thus, there is equality if and only if 

\ogf(x) = bix 2 + ax + di, \ogg(x) = b 2 x 2 + c 2 x + d 2 . (37) 

By coupling (37) with (36), we obtain that there is equality in (34) if and only if 
/ and g are gaussian densities, of variances ca and cb, respectively, for any given 
positive constant c. 

□ 
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The case r = 1 has been treated in Blachman [6], as part of his proof of the 
entropy power inequality (5). In this case 



denotes the Fisher information of the probability density /, and inequality (34) 
becomes 

(a + b) 2 I(f*g)<a 2 I(f) + b 2 I(g). 

We remark that the validity of (34) is not restricted to probability density functions. 
Indeed, it continues to hold for nonnegative functions of any given mass. 

Let us apply the result of Lemma 3 to control the sign of the right-hand side in 
formula (31). If we choose a 2 = ap/p', b 2 = f3q/q' in (34), then the coefficient of the 
term on the left-hand side of inequality (34) assumes the value 



V , Q , « r^5 PQ 



a 2 + b 2 + 2abr = a— f + (3-^ + 2 y/afi 



p q' y pq 

Let us introduce, for any given r > 1 the function 

T(a, (3) = (a + fi)(r - 1) - (a^ + 0± + I^ JK, 

\ p q V PI 



(39) 



It is clear that, as soon as for some values of a,/3 the function T(a,/3) < 0, the 
expression on the right-hand side of (31) is nonnegative, and the functional ^ u ,v(t) 
is increasing. In order to check its sign, consider that the function T is jointly convex, 
and it is such that, for any positive constant c 

F(ca, cj3) = cT(a, (3). 

Therefore, if a point (a = a , (3 = (3 ) is an extremal point, also the point (ca , c(3 ) 
is an extremal point, and T admits the half-line (3 a = a (3 of extremals. Since 



dT p fa I pq 

da p' Y (3 y p'q' 

by adding and subtracting the quantity pr/q' we obtain 

dT p p fa ppq p 

— = r +— r -l- — - r + -r 

da q' f » V P V PQ Q 



1 1 \ (\ 1 \ fa I pq p 



pr i -A-p [- + - - J -o\ — r + ~ r - 

\P Q J \P P J V P V P Q Q 
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Since 



one obtains 



1 1 _ 1 1 
p q> q pi 



(40) 



da y (3 y pV g' 



if the point (a, /3) belong to the half-line 



= — —a 
q' q 



(41) 



Same result is obtained if we impose the vanishing of the partial derivative of T with 
respect to j3. On the other hand, thanks to identity (40) 



q 1 p' 
P' q 



- + -)(r-l)-----2r 
p q J p q 



p q' 



q'rl 7 + iM 7 - - P' = 0. 



Hence, along the line (46), in view of lemma 3 the functional <& u ,v(t) is increasing 
with respect to t. Proceeding as in Section 2, namely by scaling u(x,t) and v(x,t) 
as in (17), we conclude that the functional will keep its maximum value as time goes 
to infinity, and 



lim $ (i) 

t— >oo 



f \ ( r \ 

I u(x)dx) I / v(x)dx) C(p,q,r), (42) 

JR J \JR J 



where 



C(p,q,r) =( [ {M q//p (x)^ * M p//q (x)^Y dx. 



l/r 



(43) 



Using that the convolution of Gaussian functions is a Gaussian function, M a * M@ 
M a+/ 3, we compute 



M^ p *Ml /q J dx. 



l/r 



pa q(3 (pa + q(3) 1 ^ r 
a i/ P f3i/q r Vr( pa + q p) 



1/2 



The choice a = q'/p, (3 = p'/q gives 

C(p,q,r) = (A p A q A r ,) l '\ 
where A m is defined by (2). This concludes the proof of the first part of Theorem 2. 
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The case in which 1/p+l/q — 1 + 1/r, but < p, q, r < 1 can be treated likewise. 
In this case the dual exponents p', q', r' are negative, and formula (31) takes the form 



1 d 

rdi 



a 



[ h r (x,t)dx = (a + p)(l-r) [ h r - 2 (x,t) (h x (x,t)f dx+ 

Jr L Jr 

^- [ h r - 1 (x,t)A(u 1/p ,v 1/q )(x,t) dx + p-r^L [ h r ~\x,t)B(u 1/p ,v 1/q )(x,t)dx 
\P I Jr 111 Jr 



4) 

To control the sign of the quantity into square brackets, we introduce now the 
function 



r(a, (3) = (a + (3)(l - r) - [ ay- + (3-- + 2^ U^-r 

\p\ W\ V \v\W\ 



(45) 



In this case 



^£ = l_ r + Z + I pq r 

da p' y (3 y \p'\\q'\ 



By adding and subtracting the quantity pr/q' we obtain as before 

dT fa I pq p 

r = 



da y f3 y 

if the point (a, (3) belong to the half-line 

p \p r 



\q'\ q 



-a. 



(46) 



This choice however implies that the right-hand side in (44) is non positive, and the 
functional ^ u ,v(t) decreases. This leads to the reverse Young's inequality (3). 

□ 



4 Conclusions 

In this paper we presented a new proof of the sharp form of Young's inequality for 
convolutions, as well as and its reverse form. For the sake of simplicity, this proof 
has been done in dimension n — 1. Looking at the details of the computations, it 
appears evident that the proof still holds in dimension n > 1, since the computations 
in higher dimension do not affect the constants both in formula (31) and (34), which 
are at the basis of the whole procedure. The main difference relays in the fact that 
the Gaussian functions are n-dimensional Gaussians, which lead to the additional 
presence of the exponent n in the sharp constant. Hence Theorem 2 leads to the 
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sharp inequality (1) in any dimension, without any additional (if not computational) 
difficulty Also, both Young's inequality and its reverse form are here derived by a 
unique well understandable physical principle, in the form of time monotonicity of 
a Lyapunov functional. 

Acknowledgment: The author acknowledge support by MIUR project "Opti- 
mal mass transportation, geometrical and functional inequalities with applications". 
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